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Abstract: We derive the Hawking radiation for Vaidya black hole in the tunneling picture 
from the corresponding single particle action by the use of the radial null geodesic and the 
Hamilton- Jacobi method (beyond semi-classical approximation). Both results are then 
analyzed and compared. 
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1. Introduction 

In 1974, Hawking startled the physics community by proving that black hole evaporates 
particles Q. It contradicts with classical general relativistic definition of a black hole, 
an object that nothing can escape from it [||], Q. Hawking's derivation was a quantum 
field theoretically. The semi-classical approach in modeling the black hole radiation as 
a tunneling effect was developed in the past decade and attracted many 

attentions [e.g. [|l6|]. This is mostly because of its bold physical picture and it is easier, 
mathematically. 

There are two ways to perform semi-classical analysis for a black hole radiation. The first 
is by the use of radial null geodesic method developed by Parikh and Wilczek [Q], Q. In 
this method, one has to get the expression dr/dt from the radial null geodesic condition, 
ds^ = = 0, for a metric that has the form ds^ = —a (r) dt^ + b (r) dr"^ + r'^dCl'^. Then, 
the obtained expression is used to calculate the imaginary part of the action. We then 
end up with the obtaining of the Hawking temperature Th = {k(3)~^, after equating the 
transition rate for a particle in a tunneling process T ~ exp [2ImS'] with the Boltzmann 
factor exp [—f3E]. 

The second method was developed by Padmanabhan et al. |p and it attracted some ap- 
plications (e.g. in 2-d stringy black holes |jl^). In the method, the scalar wave function 
is determined by the ansatz cj) (r, t) = exp [—iS (r, t) /h] where S (r, t) is the action for a 
single scalar particle. Inserting this ansatz into the Klein-Gordon equation in a gravi- 
tational background, one yields an equation for the action S (r, t) which can be solved 
by the Hamilton-Jacobi method. After obtaining the action, one can get the probabil- 
ity for outgoing and ingoing particles, Pout = |</'out|^ and Pin = \<pin\'^, respectively. The 
Hawking temperature can be obtained by using the 'principle of detailed balance' [5], 
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Pout = exp [—f5E\ Pin = exp [—f3E\^ since all particles must be absorbed by the black hole. 
An interesting improvement has recently been made by Banerjee and Majhi |^]. They 
expand the action for single particle in the power series of Planck constant, S {r,t) = 
'^n=o '^nK^Sn {r, t). By considering that, for all n, Sn {r, t) and Sn+i {r, t) are proportional 
to each other with the same proportionality, one can write S (r, t) = {1 + X]n=i Q^n^") •S'o i^-, t) 
By using the method of Padmanabhan et al., the action gives the correction value of the 
Hawking temperature. Even though the correction is in principle negligible, due to the 
very small of the Planck constant, one could regard this as an effect of quantum analysis 
in the semi-classical quantum gravity. Several developments have been made, including 
fermionic consideration and the relation of the method to trace anomaly ||T5| . 
In this paper, we consider a more general metric with the mass of black hole is time (t) and 
radius (r) dependent. A semiclassical approach for a dynamical black hole has also recently 
worked in [|^. As generally known, this subject is well described by the Vaidya metric. 
To get an exact {t — r) dependence of the Hawking temperature, we insist to work in using 
Schwarzschild like metric for the Vaidya space time |^] and impose the condition for a very 
slowly varying mass of the black hole, \dm/dr\ = \m'\ « 1 and \dm/dt\ = \rh\ « 1. In 
this approximation, the squared m' and m are well approximated to zero. This consid- 
eration will be useful later in showing the proportionality between each term of action's 
expansion. Our aim is to see the (t — r) dependence of Hawking temperature in the case 
of varying mass of black hole. We derive it both in the radial null geodesic and Hamilton- 
Jacobi method. For Hamilton-Jacobi method, we directly consider all of the expansion 
terms of action (beyond the method of Padmanabhan et al.) as Banerjee and Majhi had 
done for some metrics with constant black hole masses. 

The organization of our paper is as follows. In the second section, we will derive the Hawk- 
ing temperature for a metric with a varying mass by the use of radial null geodesic method. 
In the third section, the Hawking temperature is obtained by the Hamilton Jacobi method, 
by considering all terms in action expansion. In the last section, we give a conclusion for 
our work. For the rest of this paper, we use the unit dimension: Newton constant, light 
velocity in vacuum, and Boltzman constant, G = c = ks = 1- 



2. Hawking Temperature from Radial Null Geodesic Method in Vaidya 
Space-Time 

We start with the metric derived by Farley and D'Eath for Vaidya space-time 

ds^ = -( 7l - ^ + f 1 - — V rf-' + r'dn^- (2-1) 

\x [m) J \ r J \ r J 

In the above, the black hole mass m varies with time t and radius r,m = m{r,t). x {m) is 
an arbitrary function of mass and dVl^ = dO^ + sin'^ 9d(p^ is the metric of 2-sphere. Later, 
rather than adopting the law of mass evolution by Hawking |1C], [0], dm/dt = —C/rm?^ 



we choose a weaker condition, dm/dt = m = —C{m)/m'^, which is a partial derivative, 
in contrast to the Hawking law that is an exact differentiation. It will be shown that 
this choice leads to a better form of Hawking temperature. The variable C (m) in this 
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expression counts the number of particles emitted by the black hole with mass m. This 
variable will increase for decreasing mass and for early stage of evaporation, we could 
consider that this variable is a small valued-function because the evaporated mass is small 
compared to the total mass. In this line of thought, putting x{m) = —C{m)/m is still 
valid, at least at the early stage of evaporation. Our purpose in setting x (m) in this form is 
to reduce the variable contained in the metric and later on, for the case of fixed mass with 
respect to radius, to obtain time dependent Hawking temperature |11|, Th (t) = h/Snmit). 



The above considerations lead the metric {2A_) into the form 



ds^ = -(^l.^ (^1 _ + (^1 - dr^ + rW. (2.2) 

It has the general form 

ds^ = -F (r, t) + G (r, t)~^ dr"^ + r'^d^'^. (2.3) 

where in our case F (r, t) = m^^ (l — 2mr^^) and G (r, t) = 1 — 2mr^^. The general form 
( |2.3| ) will simplify our next calculation. Unless for specific purposes, we will write F (r, t) 
and G (r, t) as F and G, respectively, for the sake brevity. 

It turns out that the metric (^]^) has a coordinate singularity at rh = 2m which of course 
is time dependent. Painleve transformation which is used to remove coordinate singularity 
for a metric with time-like Killing vector, is also applicable in this analysis. Transforming 



dt^dt-\l^-j^dT (2.4) 



the metric (2^) changes into 



I - G 

ds^ = -F{r,t)dt'^ + 2FJ——dtdr + dr^ + r'^dn'^, (2.5) 

V FG 

and therefore no coordinate singularity is found. Thus, by such a Painleve transformation, 
it understandable that, in principle, a coordinate singularity in general relativity can be 
removed only by changing coordinate, without defining new physical condition or theorem. 
In the consideration that a tunneled particle moves on the path with no singularity, ( |2.5| ) 
will be useful in describing particle's dynamics. For radial null geodesic, ds^ = dil'^ = 0, 
the differentiation of radius with respect to time can be obtained from ( |2.5D as 

where +(— ) signs denote outgoing(ingoing) radial null geodesies. 

Near the horizon, we can expand the coefficient F and G by the use of Taylor expansion. 
Since F and G are {t — r) dependent, and we only need their approximation values for 
short distances from a point (horizon), we could apply the Taylor expansion at a fixed 
time. So, we can write 

F(r,t)|,~ F'(r,t)|^(r-rO + o((r-r;,)2)|^, (2.7) 
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and 



G (r, t)\^ ^ G' (r,„ t) |Jr - + O ((r - r^f) 



(21 



By the approximations ( |2.7D and ( |2.8| ) above, the dependence of radius to time in ( |2.6| ) can 
be approached by 

(r;„t) G' (r;„t) (r-r/,). (2.9) 

Now, we discuss the action of outgoing particle through the horizon. In the original work 
by Parikh and Wilczek Q, the imaginary action is written as 

rout fout Pr rout H 



ImS* = Im j prdr = Im J J dp[.dr = Im j j —^dr. (2-10) 



The above expression is due to the Hamilton equation dr/dt = dH/dpr\r where r and pr are 
canonical variables (in this case, the radial component of the radius and the momentum). 
As a reminder, the action of a tunneled particle in a potential barrier higher than the energy 
of the particle itself will be imaginary, pr = \J'lm (E — V). Different from discussions of 
several authors for a static black hole mass [e.g. Q, Q, |^], [|l2|, the outgoing 

particle's energy must be time dependent for black holes with varying mass. So, the dH' 
integration at ( |2.1[1| ) is for all values of outgoing particle's energy, say from zero to +E (t). 
By using approximation (2.9), we can perform the integration ( 2.10| ). For dr integration, 
we can perform a contour integration for upper half complex plane to avoid the coordinate 
singularity rh- The result is 

, (2.11) 

^F'(n,t)a'(r„,t) 

Since the tunneling probability is given by F ~ exp [— ^ImS] , equalizing it with the Boltz- 
mann factor exp [—(3E (t)] for a system with time dependent energy we obtain 

h^F'{rf,,t) G' {Th.t) 
Th = — . (2.12) 

As seen in expression ( 2.12| ), the Hawking temperature Th is time dependent. We will see 
it later that it is also radius dependent. Inserting the values of F' {rh,t) and G' {rh,t) with 
the condition of (m')^ ~ 0, one has 

hJ {m {rh, t)f - 2m' [rh-, t) (m (r/^, t) + 2m (r^, t)) 

Th = — 2 • (2-13) 

Svr (m {rh,t)) 

In the above, the mass m is written as m {rh,t) explicitly to remind us that the mass and 
its derivative m' {rh,t) are evaluated at rh (the radius of event horizon). Interestingly, for 
a black hole whose mass is only time dependent, not radius dependent, m' = 0, we get 

Th = (2.14) 

and thus, we recover to the time-dependent Hawking temperature. Further interesting in- 
vestigation would be to understand the black hole model with x{m) = —C{m)/m, however 
it will not be discussed in this paper. 
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3. Hawking Temperature from Hamilton-Jacobi Method in Vaidya Space- 
Time 



In this section, we would work in scalar field theory with gravitational background. We still 
work in the metric ( ^.iD and impose the condition m = —C {m)jvr? later. In this section, 
the condition [ftif' ~ and [x (m))^ ~ are very important in showing the proportionality 
between each expansion terms of action. To avoid confusion with F (r, i) and G (r, t) that 
have been used in the previous section, we rewrite metric (^.iD as 



ds^ = -/ (r, t) dt^ + g (r, t)"^ dr^ + r^dQ'^, (3.1) 
where / (r, t) = m(l — 2mr~^)/x(m) and g (r, i) = (l — 2mr~^^ . 

Massless scalar particles under the gravitational background g^^ obey the Klein-Gordon 
equation 



-a 



9m </> = 0. (3.2) 



For spherical symmetric black hole, we may reduce our attention only to (r — t) sector in 
the space-time, or in other words, we reduce to two dimensional black hole problems. In 
this consideration, equation (|3.2D under the background metric (^]^) simplifies to 



dU -Yfg{f9 + 9f) dt<P - \ [f'g + fg') drcP - fgdlcf, = 0. (3.3) 

In the above, f = f {r, t) and g = g {r,t). By the standard ansatz for scalar wave function 
(f) (r, t) = exp [—jiS (r, t)] , equation (^) leads to the equation for the action S (r, t) 



4 if 9 + m (^) (f ) - h (0) - ^ (f )> • ^'-'^ 

Now, our next step is to solve this equation. An approximation method can be applied by 
expanding the action in the order of Planck constant power, 

S{r,t) = So{r,t) + y] anffSn{r,t), (3.5) 

for n = 1,2, 3, .... The constant a„ is set to keep all the expansion terms have the action's 
dimension. Taking unit dimensions G = c = fcB = l, a„ would have the dimension of 
|-^j-2n -^^j^jgj^ refers to the mass. It is clear that this expansion would lead to a very long 
equation. Due to the very small value of the Planck constant, many authors |^], Q, |12|, 



neglect the terms for n > 1. This consideration is acceptable, and including higher 
terms is just adding correction for semi-classical derivation of Hawking temperature. By 
grouping all the terms into the same powers of h, we could write for some lowest rank as 

: fg {drSof - {dtSof = 0, (3.6) 

: ifgd?.So - i (d^So) + 2fg (drSo) (a,5i) + § {f'g + fg') dtSo 
+ i^^^dtSo - 2 (dtSo) (dtSi) = 0, 
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: zfgd^Si - idiSi + 2fg (drSo) (a.^a) + fg (drSiY + § {f'g + fg') dtSi 



+ 



(/g+g/) 
fg 



dtSi - 2 (dtSo) {dtS2) - {dtSi 



0, 



: ifgd^S2 - id!S2 + 2fg (drSo) {drS^) + +2/5 {drS2) - 2 (at5o) (9*53 

- 2 (a*52) + § + 9*52 + ^^^^dtS2 = 0, 



As obtained by Banerjee and Majhi for the metric that has a time-hke Kilhng vector, 



the metric (3J) also leads to such a relation: 

^0 : dtSo = iJTgdrSo, 



re : dtS2 = ±^/hdrS2, 



(3.7) 



To get the benefits of the conditions for black holes with slowly varying mass, (m')^ ~ 
and (m)^ ~ 0, the arbitrary function x (m) must be taken to be —m' Q. This mechanism 
means that we have used the dynamic equation for black hole mass as m' = C (m) /m. 
Unless these conditions, the set of equations ( |3.7| ) can not be obtained. Then one can obtain 
a general pattern for arbitrary n, thus for the terms of that is dtSn = ^^/JgdrSn- From 
the set of equations (^^), one can see that each 5„ (r, t) is proportional to Sq (r, t). By this 
evidence, we can write the expansion ( |3.5D into 

5 (r, t) = So (r, t) + a„/i"5o (r, t) = (l + a„/i") 5o (r, t) . (3.8) 

The extra value UnK^So (r, i) in (|3.S|) can be regarded as the correction term of the 
semi-classical analysis. So, our next step is to find the solution of Sq (r, t) satisfying dtSo = 

In the standard Hamilton-Jacobi method, Sq (r, t) can be written into two parts, the time 
part which has the form of Et and the radius part Sq (r) which is in general a radius 
dependent only. Since our metric coefficients are both radius and time dependent, the 
standard method would not be applicable. We could generalized the method by making 
an ansatz 

t 

5o(r,t) = j E{t')dt' +~SQ{r,t). (3.9) 



At the first sight, it seems that the ansatz is rather strange, that is 5 (r, t) and Sq (r, t) are 
both t and r dependent. The term J E {t') dt' is more understandable, since the emitted 
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particle's energy is continuum and time dependent. Let see how it works. 
From (13.91), one can write that 



dtSoir,t) = E{t) + dtSo{r,t) (3.10) 

and 

drSo{r,t) = drSo{r,t). (3.11) 
Since 5o (r, t) is t and r dependent, one can write 

= drSo (r, t) + dtSo (r, t) ^. (3.12) 
ar dr 



Ehminating dt/dr by the use of dr/dt = -^yfjg^ equation ( 3.12 ) can be written as 

'^^^p^ ^ ^dt~S^ (r, t) = drSo (r, t) . (3.13) 
dr ^/fg 

Combining the first equation of (|3.7D , with equations ( ^.10|) and (|3.13| ) where we should note 
that the action equation — (fg) ' dtSo (r, t) = drSo (r, t) is belong to outgoing particle 
and with the radial evolution is dr/dt = y/Jd, then we could write 

T if 9)-'^' [e it) + dtSo (r, t)) = ^^^^ T ifg)-'^' dtSo (r, t) . (3.14) 
From (|3.14| ), we can get the exact differentiation of Sq {r,t) 

Tifgy'^'Eit), (3.15) 



dSo {r,t) _ _ 1/2 



dr 

and the solution of (r, t) can be obtained by integration 



S,{r,t) = ^E{t) I (3.16) 



After inserting the mass evolution equation, m = —C {rrL)jrr?^ and the arbitrary function 
X (m) = —C {m)/m, one may identify that / and g are exactly equal to F and G stated in 
the previous section. By this equality, the integration (l3.16| ) can be evaluated by adopting 
the value of J {FG)^^^"^ dr as in obtaining expression ( |2.11[ ) from ( |2.10| ) along with it's 
approximation method (near horizon Taylor expansion). The result for the integration 
( |3l6D is 

S,ir,t) = TE{t)-^. (3.17) 
The equation gives the complete action 

S (r, t) = (l + anh-) Ue {t') dt' T E {t) -^L=\ . (3.18) 
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The signs + (— ) in expression ( |3.18| ) refer to the action for ingoing (outgoing) particle. 
Back to our first ansatz for scalar wave function, (p = exp [^^(r, t)], the wave function 
for ingoing and outgoing massless scalar particle can be read of as 



4>in (r, t) = exp 



(3.19) 



and 



'-'out 



(r, t) = exp 



(3.20) 



respectively. Consequently, from ( 3.1S| ) one can get the ingoing probability of particle as 
below 



Pin = exp 



(3.21) 



This ingoing probability must be equal to unity since all particles including the massless 
one are absorbed by the black hole. This consideration gives us the relation 



Im I E(t') dt' 



ttE (t) 



which leads the outgoing probability 

Pout = exp 



1 + > OLnK 



(3.22) 



Finally, to get the Hawking temperature from the outgoing probability ( 3.22 ), we equate 
this probability expression with exp [—(iE (t)] which in |5| is called 'detailed balance' prin- 
ciple. It yields 



'-H' 



4vr(l + En"n^") 



(3.23) 



We use the index H' to distinguish the result ( |3.23D with that of (|2.12|) . Neglecting the 
correction term X^^On^" in the denominator of ( 3.23 ), one has recovered result ( 2.12| ) 
which has been obtained by the radial null geodesies at the previous section. 



4. Conclusion 

We have worked out the Hawking temperature for the metric with no time-like Killing vec- 
tor and with the time and radius dependent coefficients. It has been shown that by the use 
of mass evolution dm/dt = —C{m)/m? and mass dependent function x{m) = —C{m)/m, 
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the resulting Hawking temperature coincides with that for the widely known time depen- 
dent one in the case of radius independent mass. By neglecting the cunK^ term in the 
denominator of the resulting Hawking temperature for the case of beyond semi-classical 
approach as given in Section 3, one recovers the result in the section 2. However, the equal 
obtained results in ( p. 12 ) and (3.22) with neglected correction terms are actually based 
on some subtleties. In section 2, we use the r that is obtained from the coordinate after 
Painleve transformation which differs with a factor half with untransformed one (directly 
obtained from Schwarzschild like metric (^.3|)). This is clear since the particle moves in the 
path with no singularity^. Thus, in section 3 we use the r in deriving the action's solution 
( |3.13 ) that is obtained from Schwarzschild like metric (without Painleve transformation), 
since the equation of action comes from it ( |3.3| ). It needs further investigations to verify 
whether the results are still the same in the suitable Painleve coordinate in deriving the 
wave equation ( p.2[ ) (parallel to section 3.2 |^]). 

The slowly varying mass with respect both to time and radius, considered in Section 3, 
affects our result on the Hawking temperature expression and the proportionality between 
each terms of action's expansion. We have restored the x{m) as —m' in section 3 to get the 
proportionality by using our slowly varying mass condition. This consideration of course 
yields an equation for mass dynamics m' = C (m) /m which needs further investigations. 
Of course, our analysis is not valid for the case of highly radiated black holes. There are 
several models that have been proposed for describing black hole radiation. But, the lack 
of experimental data does not enable one to compare between the models. 
Further investigation on the relationship between ingoing and outgoing probabilities for 
the case of black hole with varying mass/energy by the use of path integral method is 
intriguing. This work was first performed by Hartle and Hawking ||l^ for a black hole 
with constant mass and constant energy of outgoing particle. The relation proposed in the 
paper. Pout = exp [—PE (t)] Pin, is an intuitive manner as it might not be that simple. 
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